Molecular Plant

Flag leaf angle (FLA) critically determines rice yield potential under dense planting conditions. As two pivotal phytohormones determine rice FLA, the antagonistic interaction between brassinosteroid (BR) and auxin remains largely uncharacterized. We here demonstrate that BR signaling reduces auxin biosynthesis to control FLA via regulating the biosynthesis of secondary cell wall (SCW). Genetic evidence demonstrates OsYUC8 mutants disrupt SCW formation, leading to increased FLA. At the molecular level, the BR-related transcription factor OsBZR1 directly binds to and represses OsYUC8 promoter activity, thereby fine-tuning auxin-mediated SCW biosynthesis. Field evaluations reveal that osbzr1 mutants display optimized flag leaf architecture and improved yield performance under high-density cultivation. Our study not only delineates the antagonistic BR-auxin interaction governing FLA but also establishes a genetic strategy for manipulating crop architecture to maximize yield in dense planting conditions. One-sentence summaryOsBZR1-mediated repression of OsYUC8-driven auxin biosynthesis modulates secondary cell wall formation to optimize flag leaf architecture, providing a genetic strategy for maximizing yield under high-density cultivation.

Protein import across the chloroplast outer membrane is mediated by the TOC complex and regulated by GTP hydrolysis at two receptor GTPases, large and small. The former are the {beta}-barrel-containing Toc159 family of proteins (Sun and Jarvis 2023; H. Li and Chiu 2010; Shi and Theg 2013a; Kessler and Schnell 2009). The smallest member of this family in higher plants is Toc90, which has long been regarded as functionally unimportant due to the absence of obvious phenotypes in early mutant studies (Kubis et al. 2004; Infanger et al. 2011a). It was hypothesised that Toc90 is a relic of an ancestral Toc159 protein from early green lineage organisms (Hofmann and Theg 2003; Reumann et al. 2005). Here, we challenge this view by combining genetic, biochemical, phylogenetic, and in silico structural analyses to show that Toc90 is not an ancient remnant, but rather a neofunctionalised derivative of Toc159, and likely the most recently evolved member of the family. We show that Toc90 physically associates with Toc33 as a part of a photosynthetic TOC complex. Furthermore, we identify Toc90 as a regulator of chloroplast ageing. Our findings redefine Toc90 as an active and specialised component of the chloroplast protein import machinery, playing a pivotal role in senescence, especially in reducing photosynthetic protein abundance. Together, these results highlight Toc90 as a promising target for manipulating chloroplast longevity and plant senescence, with potential to improve crop productivity and yield.

Tomato (Solanum lycopersicum) is classified as a day-neutral plant, whereas its wild relatives exhibit delayed flowering under long day (LD) conditions due to higher activity of the SELF-PRUNING 5G (SP5G). In Arabidopsis thaliana, CONSTANS (CO) activates FLOWERING LOCUS T (FT), a homolog of SP5G, but whether and how CO integrates with SP5G in tomato flowering was unclear. Here, we demonstrate that SlCOL1 (the tomato CO homolog) delays flowering by directly activating SP5G in a PHYTOCHROME B1 (PHYB1)-dependent manner. Importantly, genetic and molecular analyses combining a photoperiod-responsive tomato line carrying the wild SP5G allele from S. pennellii, together with SlCOL1 and flowering-pathway mutants, revealed synergistic crosstalk among the photoperiodic SlCOL1-SP5G module, age-dependent pathway (mediated mainly by the microRNA156-SlSBP module), gibberellin (GA) pathway, and SINGLE FLOWER TRUSS (SFT) pathway. Mechanistically, we show that SP5G forms a complex with miR156-targeted SlSBP13 to directly regulate SFT expression, and that GA may interfere with SP5G activity. Together, these findings revealed a coordinated network that integrates multiple flowering signals to modulate both shared and pathway-specific targets. Our findings provide a significant advance in understanding the molecular regulation of tomato flowering and offer promising avenues for breeding strategies optimized for diverse environmental conditions and latitudes.

Plants have remarkable capacity to reconstruct entire organ systems from tissue explants. In Arabidopsis two-step tissue culture system, pluripotency regulators are specifically expressed in the middle-cell layer of the stratified callus tissue. However, regulatory mechanisms underlying the radial patterning of callus remained unclear. Here, we found that very-long-chain fatty acids (VLCFAs) synthesized in the epidermis-like outermost layer are essential for pluripotency acquisition and successful shoot regeneration. Our genetic and transcriptomic analyses revealed that the regulatory roles of VLCFAs on pluripotency acquisition involve inter-cell layer signaling in callus tissue, while they are at least partly independent of ATML1/PDF2 functions and cuticular wax synthesis in the outermost layer. VLCFAs spatially restrict procambium cell identity by non-cell-autonomously suppressing cytokinin signaling, thereby allowing for establishment of the middle-cell layer. We propose that the inhibitory relationships between layer-specific regulators underlie the intricate balance of cellular fate determination in pluripotent callus.

The Arabidopsis PWWP-DOMAIN INTERACTOR OF POLYCOMBS1 (PWO1) and Telomere Repeat-Binding Proteins 1-3 (TRB1-3, TRBs) associate with distinct and shared protein complexes involved in epigenetic regulation, yet their cooperative roles in chromatin control and plant development remain largely unexplored. Here, we show that the interaction between PWO1 and TRBs is evolutionarily conserved. Both PWO1 and TRBs associate with plant telomeres, interact at these regions, and are co-enriched at subsets of interspersed telo-box motifs across regulatory regions genome-wide. TRBs facilitate PWO1 binding at shared genomic regions, including telo-box motifs. PWO1 and TRBs share a substantial number of genomic targets and preferentially bind chromatin regions associated with transcriptionally active states, whereas TRBs alone associate with repressive marks at thousands of loci. Genetic analyses show that the pwo1 trb1 trb3 triple mutant displays severe developmental defects, including main stem arrest and early maturation associated with aberrant lignin deposition in interfascicular tissues. In the triple mutant, key enzymes in the lignin biosynthesis pathway are upregulated, indicating that PWO1, TRB1, and TRB3 cooperatively regulate secondary cell wall formation. Together, our findings provide new insights into how PWO1 and TRBs cooperate to regulate chromatin states and orchestrate plant development, highlighting their central role in controlling gene expression programs. Significance statementThis study shows that PWO1 and TRB proteins co-occupy telomeres, including interspersed telo-box motifs, to regulate chromatin organization and plant development, particularly ectopic lignin deposition. Our findings reveal how these nuclear protein factors coordinate epigenetic states in Arabidopsis thaliana, providing a framework for understanding the control of developmental programs. O_FIG O_LINKSMALLFIG WIDTH=198 HEIGHT=200 SRC="FIGDIR/small/740627v1_ufig1.gif" ALT="Figure 1"> View larger version (38K): org.highwire.dtl.DTLVardef@1370171org.highwire.dtl.DTLVardef@3faea5org.highwire.dtl.DTLVardef@e30921org.highwire.dtl.DTLVardef@16c948e_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical abstract:C_FLOATNO Evolutionarily conserved PWO-TRB interactions and their shared roles in chromatin regulation and plant development. Created with BioRender.com. C_FIG

Nitrogen (N) is a major determinant of plant growth and productivity. Because soil N availability and internal N demand fluctuate, plants have evolved sophisticated mechanisms to coordinate N acquisition and utilization at the whole-plant level. However, how this coordination is achieved remains poorly understood. Here, we show that N-inducible LATERAL ORGAN BOUNDARIES DOMAIN transcription factors LBD37, LBD38, and LBD39 (LBDs) function as repressors of local N uptake and assimilation and systemic N-demand signaling in Arabidopsis. Triple mutants lacking these three LBDs displayed enhanced nitrate influx and increased accumulation of nitrate, amino acids, and total N. Transcriptome analysis identified an array of N-starvation- and nitrate-inducible genes derepressed in shoots and roots, including C-TERMINALLY ENCODED PEPTIDE (CEP) and CEP DOWNSTREAM (CEPD) genes, as well as genes involved in N uptake and assimilation. Grafting and genetic analyses revealed that LBDs gate the systemic N-demand signaling relay by repressing CEP and CEPD expression organ-autonomously. We also found that LBDs locally repress genes involved in N uptake and assimilation through a distinct regulatory mechanism. We propose that LBDs are key transcriptional repressors in a regulatory framework for optimizing N acquisition and utilization under fluctuating N conditions at the whole-plant level.

The ATP-binding cassette (ABC) transporter ABCB1 transports both auxins and the brassinosteroid brassinolide (BL). ABCB1-mediated auxin (IAA) transport depends on interaction with the FKBP42 protein TWISTED DWARF1 (TWD1), yet the mechanism underlying substrate selectivity remains unclear. Here, we confirm dual IAA and BL transport by ABCB1 and show that the two substrates compete for ABCB1-mediated transport. We demonstrate that a conserved proline residue (P1008) is required for IAA, but not BL, transport. Furthermore, we identify TWD1 as a calmodulin-activated peptidyl-prolyl cis-trans isomerase (PPIase) that selectively enhances IAA, but not BL, transport through isomerization of the E1007-P1008 peptide bond. Loss of TWD1 PPIase activity abolishes ABCB1-mediated IAA export without affecting BL transport, revealing a regulatory mechanism that prioritizes auxin over brassinosteroid transport. Cryo-EM, molecular docking, and molecular dynamics simulations of wild-type ABCB1 and the ABCB1P1008G variant suggest that the P1008 loop mediates long-range communication between the nucleotide-binding domain surface and the substrate-binding pocket, thereby providing a structural mechanism by which TWD1 modulates substrate specificity. In summary, our findings establish the structural basis of ABCB1 substrate selectivity and uncover a post-translational mechanism that selectively regulates auxin transport while preserving brassinosteroid transport. One Sentence SummaryTWD1 prioritizes ABCB1-mediated auxin over brassinosteroid transport

Protein-protein interactions (PPIs) play a crucial role in enabling proteins to carry out their functions within various biological processes (Hui et al., 2003). Since the introduction of the yeast two-hybrid (Y2H) method for PPI detection in 1989 (Fields and Song, 1989), the identification of PPIs has become a significant focus in modern biological research. PPI goes beyond examining individual proteins, allowing researchers to establish a comprehensive network that regulates biological processes. Rice, as a key model organism in plant biological studies, has been at the forefront of PPI research. In 2008, prominent rice scientists in China called for concerted efforts to define a comprehensive protein-protein interaction network experimentally, which aimed to facilitate the prediction of the functional mechanisms operating throughout a plants lifecycle (Zhang et al., 2008). With efforts for 2 decades, the experimentally identified rice PPIs have reached over ten thousand. Several public databases have been established to systematically collate and store PPIs, including STRING (Szklarczyk et al., 2019), BioGRID (Oughtred et al., 2020), IntAct (del Toro et al., 2022), PRIN (Gu et al., 2011), RicePPINet (Liu et al., 2017) and RiceNet v2 (Lee et al., 2015). However, most PPI datasets in rice stem from computational predictions, while experiment-based rice PPI datasets are fragmented due to the lack of systematic profiling at the rice PPIome level, which largely hinders information sharing in the rice research community. To bridge this gap, we constructed the Port of Protein-Protein Interactomes (POPPIN; https://riceome.hzau.edu.cn/poppin/), an integrated database dedicated to sharing experimentally verified PPIs and functional clues in rice. Empowered by high-throughput PPIome profiling technologies and text mining assisted by a large language model (Huang et al., 2025; Liu et al., 2025), POPPIN currently has deposited over 150,451 pieces of rice PPI-related information. Additionally, POPPIN provides detailed protein information, including GO annotations, subcellular localizations, domains, trait ontology (TO) information, and hyperlinks to external biological databases. Through offering a user-friendly web interface for search and dynamic network visualization, POPPIN serves as the first large-scale, experiment-based database for searchable PPIs in rice, and has the potential to be extended to other species under this structural framework.

Under nitrogen limiting conditions, legume plants interact with nitrogen fixing bacteria known as rhizobia, resulting in the formation of a new organ, the nodule. This process is accompanied by dramatic changes in gene expression, which operate at different levels. A previous study revealed that histone methylation is differentially modulated during nodulation. However, the histone methyl transferases and demethylases involved in this modulation have not been characterized. In this study we report the identification of the Medicago truncatula putative histone lysine demethylase MtPKDM9B, which is subject to alternative splicing (AS), and the differential modulation of AS variants at translational level during nodule symbiosis. Knockdown of MtPKDM9B impaired infection by rhizobia, nodule development, bacterial viability and the expression of the leghemoglobin coding gene MtLHB1. MtPKDM9B is the putative ortholog of Arabidopsis EARLY FLOWERING 6 (ELF6/AtPKDM9B) gene involved in the removal of the repressive mark H2K27me3. A combination of ChIP-seq and RNA-seq experiments revealed that MtPKDM9B is required for demethylation of H3K27me3 in regions nearby or contained within gene bodies of symbiotic genes and the upregulation of the cognate mRNAs in response to rhizobia, including those encoding the putative ubiquitin ligase MtPUB2, the MYB transcription factor MtMYB040 and the auxin conjugating enzyme MtGH3 (Gretchen Hagen 3). Our findings illustrate how AS and translational regulation of this plant specific histone lysine demethylase contributes to the removal of the repressive mark H3K27me3, promoting transcriptional activation of symbiotic genes required for the formation of functional nitrogen fixing nodules.

Coordinated organogenesis and rhythmic stem development determine shoot morphogenesis, which governs distinct plant architecture and crop yield. The shoot apical meristem (SAM) drives organ initiation and stem growth, with its activity sustained by the spatially defined CLAVATA3 (CLV3)-WUSCHEL (WUS) feedback loop. However, how the temporal CLV3-WUS module modulates SAM activity remains elusive. Here, we demonstrate that both CLV3 and WUS are required for circadian stem elongation. The clock component TIMING OF CAB EXPRESSION1 (TOC1) directly activates and sustains the circadian expression of CLV3, and this regulatory process is controlled by the APETALA1-CIRCADIAN CLOCK-ASSOCIATED1 module. Notably, the rhythmic expression pattern of CLV3, rather than its transcript abundance, plays a more important role in regulating WUS expression and sustaining SAM activity, and CLV3 confers circadian rhythmicity to WUS protein accumulation. Transcriptome and ChIP seq analyses reveal that WUS modulates multiple developmental programs, including meristem homeostasis, phytohormone responses, and cell wall remodeling. Mechanistically, WUS directly activates PECTIN METHYLESTERASE5 expression to promote stem elongation. Collectively, our findings uncover a spatiotemporal regulatory cascade that fine-tunes SAM activity to coordinate rhythmic organogenesis and stem development, thereby establishing optimal phyllotactic patterning in plants.

Plant regeneration underpins plant developmental plasticity, tissue culture, genetic transformation and crop improvement. Although numerous related omics datasets have been generated, specialized multi-omics databases for this field are still scarce. Here, we constructed PlantRegMoD, an AI-powered integrated database dedicated to plant regeneration. This platform hosts 20.54 TB standardized multi-omics data from 147 projects across 32 plant species and 2,593 samples. We established a unified hierarchical classification system covering five major categories and nine regeneration models, and curated 236 regeneration genes as well as their 28,190 homologs across 58 representative plant species. It contains extensive transcriptomic resources across all regeneration models, together with 196,423 single cells and over 8.81 million epigenetic peaks to dissect cellular heterogeneity and multi-layered epigenetic regulation. Equipped with nine online omics-related tools and a RAG-based intelligent Q&A system, PlantRegMoD greatly reduces bioinformatic barriers and serves as a robust resource for mechanistic, functional and evolutionary studies of plant regeneration.

Global warming is increasing the frequency and intensity of high-temperature episodes, limiting plant productivity. However, the molecular mechanisms integrating primary metabolism with the heat stress response remains poorly understood. Here, we show that lysine biosynthesis contributes to the coordination of physiological, metabolic and transcriptional responses to heat stress in Arabidopsis thaliana. We compared wild-type, the lysine-biosynthesis mutant dapat, and the salicylic acid (SA)-biosynthesis and signaling mutants sid2-1 and npr1-3 under prolonged warming (6C above control for 7 days) and heat shock (38C for 6 h), followed by recovery. We assessed growth, gas exchange, photosynthetic performance, free SA, salicylic acid glucoside (SAG), salicylic acid glucose ester (SGE), and total SA content, primary metabolite profiles, heat-stress-responsive gene expression and transcriptome-wide changes by RNA sequencing. Before heat stress, dapat mutant presented a distinct metabolic state, marked by amino-acid accumulation, altered organic-acid profiles, reduced soluble sugars and elevated endogenous SA. This metabolic configuration persisted during prolonged warming, whereas WT and SA-pathway mutants underwent more dynamic reprogramming. Heat shock, by contrast, elicited a more convergent response across genotypes. Despite reduced basal PSII efficiency, dapat maintained photosynthetic performance during prolonged warming and recovered. Its transcriptional response, however, differed from that of WT and SA-pathway mutants: selected heat-responsive genes were constitutively or more strongly expressed, whereas some canonical heat-stress regulators showed weaker induction after heat shock. RNA-seq further revealed a largely conserved core heat-shock response but genotype-dependent regulation of defense, hormone and amino-acid-metabolism programs, particularly during recovery. Together, these findings indicate that impaired DAPAT activity establishes a metabolically primed but energetically constrained state that reshapes gas exchange, photosynthetic acclimation and heat-responsive transcription. Lysine homeostasis therefore emerges as a regulatory node linking primary metabolism and SA accumulation with SA-dependent and SA-independent components of heat-stress acclimation.

Brown spot (BS), caused by the fungal pathogen Bipolaris oryzae, is a major disease threatening global rice production. However, the genetic basis of host BS resistance remains unclear. Here, we identified brown spot resistance 1 (bsr1), a quantitative trait locus conferring BS resistance, by map-based cloning. We show that bsr1 encodes a sucrose transporter and that a near-isogenic line carrying bsr1 (bsr1-NIL) in the susceptible Koshihikari genetic background exhibited resistance to BS by suppressing sucrose efflux into the apoplast after pathogen attack. Furthermore, bsr1-NIL also showed strain-specific resistance to bacterial blight caused by Xanthomonas oryzae pv. oryzae through the same mechanism. These findings demonstrate that bsr1 confers dual resistance to fungal and bacterial diseases by regulating sucrose efflux. Our study identifies a previously unrecognized mechanism underlying resistance to both BS and bacterial blight and highlights bsr1 as a promising target for breeding disease-resistance rice cultivars. Rice (Oryza sativa L.) is a staple food for more than half of the worlds population1. Brown spot (BS), caused by the fungus Bipolaris oryzae, is one of the most prevalent fungal diseases of rice, and its incidence has increased under global warming2. BS infects coleoptiles, leaves, leaf sheaths, panicle branches, glumes, and spikelets, and severe infection can substantially reduce grain yield.

Groundnut (Arachis hypogaea) is an important oilseed crop, and its production is severely affected by drought worldwide. Improving the drought-tolerance of groundnut without compromising its agronomic traits is a major challenge for researchers. Here, a meta-transcriptome analysis identified tandemly duplicated G{beta} subunits (AhG{beta}1 and AhG{beta}3) as early signaling genes under water deficit conditions. Heterotrimeric G-protein {beta}-subunit is a major signaling element that can regulate diverse responses in plants. Heterologous complementation studies in Arabidopsis G{beta} mutant (agb1-2) demonstrated evolutionary conservation of AhG{beta} function, while the transgenic overexpression and RNAi knockdown in groundnut established AhG{beta}1 as a positive regulator of plant growth, yield, and drought tolerance. Integrated transcriptomic, transcription factor regulatory network, and metabolomic analyses demonstrated that AhG{beta}1 coordinates drought-responsive transcriptional and metabolic reprogramming characterized by suppression of cell wall-associated processes and activation of phenylpropanoid, flavonoid, and isoflavonoid biosynthesis. This reprogramming resulted in enhanced accumulation of antioxidant flavonoids, including quercetin and rutin, reduced ROS accumulation, and improved drought tolerance. Furthermore, exogenous application of these flavonoids enhanced drought recovery, highlighting their role as protective downstream metabolites. Collectively, our findings establish AhG{beta}1 as a conserved regulator linking early G-protein signaling with transcriptional and metabolic reprogramming to maintain redox homeostasis during drought stress and identify promising molecular and metabolic targets for improving drought resilience in groundnut.

Fluctuating light strongly limits photosynthetic efficiency in plants, yet the mechanisms coordinating stomatal dynamics with rapid photosynthetic induction remain incompletely understood. Here, through a forward physiological screen in Arabidopsis, we identify the abscisic acid (ABA) transporters NPF2.4 and NPF2.5 as regulators of photosynthetic responses to fluctuating light. Loss-of-function npf2.4 and npf2.5 mutants exhibited accelerated stomatal opening and photosynthetic induction during light transitions, resulting in enhanced biomass accumulation and seed yield without compromising long-term drought tolerance. NPF2.4 and NPF2.5 mediated ABA uptake in yeast and were predominantly expressed in the root pericycle and cortex/epidermal tissues, respectively. Reciprocal grafting, xylem sap ABA measurements, and ABA treatment to the roots demonstrated that both transporters regulate root-to-shoot ABA transport. Our findings reveal that root-derived ABA modulates stomatal dynamics and photosynthetic efficiency under fluctuating light, establishing long-distance hormone transport as a key regulator of plant carbon gain in dynamic environments.

Topologically associating domains (TADs) are essential structural and functional modules of the genome that play a crucial role in regulating gene expression. In this study, we systematically investigated the conservation and evolution of TADs in five closely related crops, including maize, sorghum, coix, foxtail millet and broomcorn millet. Our results show that 74% of TAD boundaries are conserved between two inbred maize lines, B73 and Mo17, and that approximately 50% or more of TAD boundaries are conserved across different crop species. TAD number remains relatively stable in the face of changes in genome size. However, the length of TADs varies depending on genome size. Furthermore, we found that large-scale transposable element expansion leads to TAD expansion, while chromosomal inversions lead to TAD fusion and the formation of new TAD boundaries. Frequent chromatin interactions between subgenome chromosomes occur after whole-genome duplication. Moreover, we also found that crossovers are enriched at TAD boundaries in maize, indicating the importance of TADs as a fundamental unit during species evolution. Overall, our study provides insights into the conservation and evolution of TADs in crop genomes and their roles in genome organization and function.

Iron (Fe) deficiency is a major constraint for plant growth and triggers extensive physiological and transcriptional reprogramming to maintain Fe homeostasis. Here, we identify the glucosinolate-associated transcription factors MYB28 and MYB29 as previously unrecognized regulators of Arabidopsis thaliana adaptation to Fe deficiency. Across various growth systems, loss of MYB28 increased sensitivity to Fe deficiency, whereas the myb28myb29 double mutant displayed stronger chlorosis, reduced root growth and impaired biomass accumulation, indicating cooperative but unequal functions of these transcription factors. Despite their enhanced Fe-deficiency phenotype, double mutant plants accumulated higher Fe levels in roots and exhibited stronger induction of canonical Fe-deficiency responses, suggesting impaired Fe utilization or distribution rather than defective Fe uptake. RNA-seq revealed extensive transcriptional reprogramming under Fe deficiency, with pronounced deregulation of genes involved in Fe homeostasis, redox processes and growth, particularly in the double mutant. Among these, SCOPOLETIN 8-HYDROXYLASE (S8H) emerged as a major target gene of MYB28. The expression of S8H was almost abolished in myb28 and myb28myb29 mutants, whereas expression of other coumarin biosynthetic genes remained largely unaffected. Promoter activation assays demonstrated that MYB28 activates the S8H promoter, and metabolic analyses showed accumulation of scopolin together with reduced fraxin levels in the mutants, consistent with impaired S8H activity. Collectively, our results identify MYB28 as a key regulator linking specialized metabolism to Fe homeostasis through control of coumarin biosynthesis, thereby expanding the biological functions of MYB28 and MYB29 transcription factors.

Manilkara zapota (Chiku) is a tropical fruit tree in the Sapotaceae family. In this study, we assembled the first Telomere-to-telomere (T2T) genome assembly of M. zapota, spanning 1.6 Gb with high contiguity (N50 = 134 Mb). Comparative analysis revealed a very recent divergence and a shared WGD among Sapotaceae members. A newly constructed Ancestral Sapotaceae Karyotype (ASK) suggests a fusion-dominated chromosome structure in the family. The genome exhibits lineage-specific proliferation of transposable elements (TEs) accompanied by a recent burst of LTR-Gypsy elements (~0.09 mya). Genome-wide methylation profiling revealed elevated CpG and CHG methylation levels and age-dependent reinforcement of TE silencing, suggesting strong epigenetic regulation that stabilises the repeat-rich genome while maintaining gene expression. Collectively, these findings advance our understanding of repeat-rich genome evolution in plants and provide a valuable framework for genomic-assisted breeding in M. zapota.

Climate change exposes plants to recurrent and interacting stresses, yet the extent to which these effects persist across generations, and the mechanisms involved, remain unclear. We propagated Arabidopsis thaliana wild type (WT, Col-0) and nitric oxide homeostasis mutant gsnor1-3 (hereafter, gsnor-ko) for five successive generations. Plants were grown under control, drought, elevated CO2, O3, warm temperature, and combined treatment scenarios for the first three generations (G1-G3), followed by two recovery generations under control conditions (G4-G5). We quantified rosette growth, photosynthetic traits, seed production, and transcriptome dynamics by RNA-seq. Across environments, gsnor-ko showed reduced vegetative growth and reproductive output relative to WT. Transcriptomic responses were strongly scenario- and generation-dependent, with the largest differential expression shifts observed under warm-climate and combined-treatment conditions. Compared with WT, gsnor-ko displayed broader gene overlap across generations and stronger retention or reconfiguration of stress-responsive states after stress withdrawal. Functional enrichment and candidate-gene analyses identified pathways/components linked to DNA methylation, heterochromatin maintenance, histone ubiquitination, m6A RNA regulation, and methyl-donor metabolism. Together, these results support a model in which GSNOR activity promotes transcriptomic recovery after repeated climate stress, whereas impaired GSNOR function shifts responses toward multi-generational persistence and epigenetically associated regulatory reconfiguration. HighlightGSNOR-dependent nitric oxide homeostasis promotes transcriptomic resetting after repeated climate stress, whereas impaired NO homeostasis favours multigenerational persistence and chromatin- and RNA-linked regulatory reconfiguration. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=118 SRC="FIGDIR/small/742973v1_ufig1.gif" ALT="Figure 1"> View larger version (31K): org.highwire.dtl.DTLVardef@15aa6c2org.highwire.dtl.DTLVardef@7a419eorg.highwire.dtl.DTLVardef@69f278org.highwire.dtl.DTLVardef@1adbc16_HPS_FORMAT_FIGEXP M_FIG C_FIG

Precise quantitative and spatial control of plant gene expression is typically attributed to cis-regulatory elements in promoters and other non-coding sequences. The extent to which exon-encoded regulatory information shapes plant gene expression and how remains unclear. Here we used BRASSINOSTEROID INSENSITIVE 1 (BRI1) as a model locus and found that it contains an exonic enhancer (EE), which is essential to establish the spatial BRI1 expression pattern. Its disruption restricts BRI1 expression and impairs BRI1-dependent root growth, linking exon-encoded regulatory information and development. The BRI1 EE acts in concert with its promoter, and related receptor kinase genes likewise harbor regulatory information within their coding sequences. Genome-wide surveys show that candidate EEs are associated with highly connected chromatin interaction networks and transcriptionally active loci. Moreover, species comparisons suggest that common features of exon-encoded regulatory information extend across diverse plant lineages. Thus, exon-encoded regulatory information is a vital component of plant gene regulation.